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ABSTRACT 

The  system  studied  is  a  surface-to-air  guided  missile  whose 
aerodynamic  and  control  characteristics  can  be  represented  by  a  set 
of  linear  differential  equations.     The  objects  of  the  study  are  to  select 
a  set  of  control  system  gains  which  will  give  specified  performance  of 
the  missile  within  the  roll  and  yaw  subsystems  and  at  the  same  time 
to  minimize  the  coupling  effects  whereby  motion  about  one  of  the  axes 
causes  dynamic  response  about  the  other  axis. 

The  system  is  studied  analytically  using  the  mathematics  of 
control  system  synthesis  and  design.     Analytic  investigations  are 
paralleled  by  system  simulation  on  a  REAC  analogue  computer. 

Using  the  above  methods,   gains  are  selected  which  give 
specified  response  in  the  roll  stabilization  and  yaw  control  subsystems 
and  suitably  damped  response  due  to  coupled  motions  between  systems. 
The  study  also  investigates  the  possibilities  of  compensation  to 
attenuate  or  eliminate  coupled  response. 
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Title:  Associate  Professor  of 
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IV 


OBJECT 


The  object  of  this  paper  is  to  analyze  a  surface-to-air  missile 
using  equations  which  describe  its  aerodynamic  and  control  system 
responses;  to  select  values  for  control  system  parameters  which  will 
give  specified  response;  and  to  investigate  the  possibility  of  compen1 
sation  which  will  eliminate  coupled  response  effects  between  sub- 
systems. 


CHAPTER  1 


INTRODUCTION 


1. 1    Physical  Description  of  the  Missile  System 

The  missile  studied  is  a  surface  launched  antiaircraft  weapon. 
It  is  a  homing  missile  using  proportional  navigation.     Thus  the  missile 
borne  receiving  antenna  tracks  the  target  and  the  guidance  system 
continually  measures  the  angular  rate  of  change  of  the  line  of  sight  to 
the  target.     The  guidance  computer  generates  command  signals  for 
yaw  and  pitch  autopilot  systems.     The  autopilot  systems  react  by 
changing  the  missile  heading  at  a  rate  proportional  to  the  rate  of  change 
of  the  line  of  sight.     Because  of  cruciform  missile  shape,    the  yaw  and 
pitch  autopilots  have  the  same  form.     Polarization  of  the  receiving 
antenna  and  the  insertion  of  additional  autopilot  commands  in  earth 
coordinates  require  roll  stabilization  of  the  missile. 

The  general  control  problems  of  the  autopilot,  in  response  to 
commands  generated  by  the  guidance  computer,  may  be  stated  as: 

1.  An  output  acceleration  in  yaw  equal  to  the  command 
acceleration  in  yaw  with  as  fast  a  response  as  possible. 

2.  Same  response  as  above,  but  in  pitch. 

3.  Maintenance  of  roll  stabilization  (zero  roll  angle)  in  the 
presence  of  induced  roll  torques. 

1.  2    System  Equations 

The  coordinate  system  to  be  used  is  shown  in  Figure  1.1. 
The  axes  are  centered  at  e.g.   and  rotate  with  the  missile  airframe. 


Fig.   1.  1    Missile  Axis  Orientation 


Consider  i,  j,   k  to  be  unit  vectors  along  the  xj,  z  axes  re- 


spectively.    Define  the  following: 
F       =        applied  force     = 


V 
Z3* 

m 
M 
H 


'1  +  Ff+F    $ 


X 


velocity  of  the  missile    =      u  i  +  vj  +  wk: 
angular  velocity  of  the  coordinate  system 

A    J         A    ,       A 

pi  +    qj  +  rk 

mass  of  missile 

applied  moment    =    M    i  +  M    i+Mk 

rr  X  y  J  Z 

angular  momentum  of  missile 

I    pi  +  I    qj  +  I    rk 
x  ^         y  ^J        z 


(1.1) 
(1.2) 

(1.3) 
(1.4) 
(1.5) 

(1.6) 


Consider  the  missile  to  be  a  rigid  body,  and  the  axes  chosen  to 
be  principle  axes.     The  force  and  moment  component  equations  of 
motion  are: 
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M 
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m  (u  -  rv  +  qw) 
m  (v  -  pw  +  ru) 
m  (w  -  qu  +  pv) 
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If  we  consider  the  gravity  components  of  force  as  inputs  to  the 
system. which  are  trimmed  out  in  steady  flight  by  command  bias  signals 
we  may  write  the  applied  forces  and  moments  as: 


x 


where: 


M3 

M 
M 
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C 
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T  -  QSCdQ(M) 
QSCyfl3,  i     ,   6,   M) 
QSCN(<*,   i     ,   6,   M) 
QSdC^  (p,  a,  /3,    6,   i 
QSdCm  (p,  a,   6,   M) 


y 


p 


M) 


do 


QSdCn  (/3,  i     ,    6,   M) 


Rocket  thrust,   aligned  with  x  axis 

Dynamic  pressure  =ypv 

Mach  number 

Reference  area,   (usually  body  cross  section) 

Nondimensional  drag  coefficient 

Nondimensional  yawi  force  coefficient 


(1.13) 

(1.14) 
(1.15) 
(1.16) 

(1.17) 
(1.18) 

(1.19) 

(1.20) 
(1.21) 
(1.22) 
(1.23) 
(1.24) 


CN    =        Nondimensional  pitch  force  coefficient  (1.  25) 

C,      e        Nondimensional  rolling  moment  coefficient  (1.  26) 

C       =        Nondimensional  pitch  moment  coefficient  (1.  27) 

m  r 

C        =        Nondimensional  rolling  moment  coefficient  (1.  28) 

d        =        Reference  moment  arm  (usually  missile 

body  diameter)  (1.  29) 

j8        =        Yaw  angle  of  attack       =     arc  sin  rr  (1.  30) 

w 
a        =        Pitch  angle  of  attack     =     arc  sin  fr  (1.31) 

i         =        7      =      Control  surface  deflection  to  produce 
y      " 

yaw  forces  and  moments  (1.  32) 

i         a        Control  surface  deflection  to  produce  pitch 

forces  and  moments  (1.  33) 

6        =        Control  surface  deflection  to  produce  rolling 

moment  (1.  34) 

From  Figure  1.  2  we  see  that  y  ,   i      and  6  can  be  related  to  the 


actual  control  surface  deflections  fi^-v  by 


7      = 

6®    +      6® 

i      = 
P 

*©   +    6<D 

6     = 

6©    "    6© 

(1.35) 
(1.36) 

;®     +     6®  (L37) 
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Fig.   1.  2    Positive  Deflection  of  Control  Surfaces 


Now  consider  the  missile  to  be  in  the  following  steady  state 
condition: 

1.       Constant  speed  (constant  Mach)  where 

I      2~~  2~~  2  I      2   "    ...   2 

V  U         +     V         +     W  =     V  U         + 


V  =   V     = 

o  o  o  o 


io-  +   yQ~  (1.  38) 

2.       Missile  trimmed  to  constant  pitch  angle  of  attack,   a     . 


This  means  that 


M      =   QSdC       =    0 
y  m 


(1.39) 


and 


F       =  QSCtvt   =     constant 

Z  IN 


The  above  may  be  expressed  by: 


v     =     p      =r      =8      =6       =7       =0 
o         *o         o         o         o  'o 


(1.40) 


(1.  41) 


and 

i        =       constant  (1.  42) 

P 

q       =       constant  (1.  43) 

Allow  motions  of  v,   r  and  p  as  small  perturbations  about  the 
steady  state  values  noted  above,    and  consider  their  products  to  be 
negligible.     We  find,   on  substituting  into  the  equations  in  which  v,   r 
and  p  occur: 

ac  ac  ac 

=    m   [  V     (cos   8)  8  -  V     (sin  a  )  p  +  V     (cos  a  )  r]         (1.  44) 

L        O  f'f  0  0  if  0  o 

9Ci        ac,        ac,         ac,  /,   *r\ 

Mx  ■  QSd  ^  P  +  5p    *  +  B«   6   +  ^    ^  *  JxP  ' 

ac  ac  ac 

My  -  QSd  (  -f  fi  +    yJ  T    +    — n  6  )  -   Izr  -  (Ix  -  Iy)pqo  (1.  46) 

The  partial  derivatives  above  are  stability  derivatives  '     and 
are  obtained  empirically  for  airframe  and  conditions  specified  by  wind 
tunnel  tests. 

Substitute  the  below  assumptions  and  definitions  into  equations 
(1.44)  -  (1.46)  to  arrive  at  the  final  linearized  aerodynamic  performance 
equations  of  the  missile,   equations  (1.64)  -  (1.66). 

cos  8=1,   since  8  less  than  15  (1.47) 

Q  =    constant  since  V  is  constant  and   p   is  constant 

at  steady  altitude  (1.  48) 

m  ,   I     and  I      are  constants  since  propellant  burning 

is  slow  process  (1.  49) 

p     *     \     =    roll  rate  (1.  50) 

q     =     6     =     pitch  rate  (1.  51) 


r    i    tjj     •    yaw  rate  (1.  52) 

(1.53) 

(1.54) 

(1.55) 
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P  5  QSd^3  p  •  (1>57) 

x 
~  QSd  57.3    9C1  (1.58) 
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(1.59) 


L    ,     ^^    -1  (1.60) 

x  ™ 

M  -■  ^zr1  Sy  (1-61> 

o 

N    f     QSd  57.  3    ^n  (1.62) 

Iz  36 

Table  A.  1,  Appendix  A,   tabulates  the  above  aerodynamic 
coefficients  at  three  pitch  angles: 

aQ=    0°  ,  12°,   24°  (1.63) 


Yaw  force  equation: 

(A  j3  +  By  +  M5)  =  j3  -  <f>   sin  a     +  ^   cos  <*  (1.  64) 


Yaw  moment  equation: 


^r  Negligible 


(1.65) 


Roll  moment  equation: 

G6  +  F</>+  Hj3  +  L7  =  V  (1.66) 


The  above  development  is  based  on  information  received  from 
Convair  (Pomona)  and  roughly  parallels  a  more  rigorous  development 
of  linearized  aircraft  performance  equations  outlined  in  Reference  1.1, 
Sections  II  w-  I  through  II  -  8. 

Having  the  aerodynamic  performance  equations  we  need  now  to 
consider  the  control  systems  in  yaw  and  roll  (pitch  being  considered 
constant  in  this  analysis). 

The  following  sensors  are  used  in  the  missile: 

F  F 

1.  Accelerometers  which  in  effect  measure  —   and    — ">   , 

m  m 

or  with  suitable  scale  factors,   (Aj3  +  B7  +  M6)  and  (Aa  +  B7  +  M6). 

•       •       • 

2.  Rate  gyros  which  measure  $  ,   ^  ,  6  . 

The  accelerometers  provide  information  for  feedback 
comparison  with  command  signals  and  the  angular  rate  signals  are 
used  for  damping. 

The  sensor  outputs  used  in  this  analysis  are  related  to  the 
sensed  quantities  as  shown  below: 

(A3  +  Bt  +  M6)Ms)  =   G1(s)(A/3  +Bt  +  M6)(s)  (1.67) 

5>'(s)    =  G2(s)  lis)  (1.68) 

J'(s)    =  G3(s)  J(s)  (1.69) 


where  the  primed  quantities  are  the  sensor  outputs  and  the  sensors  are 
quadratic  of  the  form 

2 

CO 

G(s)  =   -2 g  (1-70) 

s     +  2;ion  s  +  un 

The  control  system  equations  for  the  autopilots  considered  in 
this  analysis  are  the  yaw  control  equation 

7c(s)=   [  Kx  (Aj3  +    BT+  M6)'(s)+ (K2+K3s)  ^'(s)]      (1.71) 

and  the  roll  control  equation 

6c(s)  =  -  [  -^  +  K5   +  K6  s]J»  (s)  (1.  72) 

These  are  integrated  through  rate  servos  to  control  surface 
deflections 

G  (s)  . 

t(s)=  -4—  ~^<s>  <L73> 

and  „ 

6(s)=    -~-6  (s)  (1.74) 

where  G.  and  Gj-  are  of  the  form  indicated  in  equation  (1.  70). 

A  functional  diagram  representation  of  the  system  is  shown  in 
Figure  1.  3. 


1.  3    Performance  Specifications 

Consideration  of  the  equations  representing  the  aerodynamic 
response  functions  of  Figure  1.  3  shows  that  interaction  exists  between 
(A/3  +  B7  +  M6),J   ,  J    and  y  and  6  at  a    >    0°. 

Since  the  system  can  be  decoupled  at  a    =  0   ,   specifications  for 
response  at  that  pitch  angle  are  to  be  met  first.     As  set  forth  by 
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Gc 
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Fig.   1.  3    Missile  System  Functional  Diagram 

Convair  in  stating  the  problem,   the  responses  to  be  achieved  (in 
uncoupled  mode)  are: 

1.  Time  constant  of  the  yaw  subsystem  to  step  command 
(A/3  +  B7)     to  be  less  than  .6  seconds.     (Time  constant  here  is  defined 
as  the  time  for  the  quantity  to  reach  6  3%  of  its  steady  state  value. ) 

2.  Time  constant  of  the  roll  subsystem  response  to  a  step 
torque  (6T)  to  be  less  than  .  5  second.     (The  time  for  <f>   to  decay  to 
less  than  37%  of  its  peak  value). 

The  only  requirement  in  coupled  modes  (a     -   12    t   24   )  is  that 
the  system  be  stable  and  reasonably  damped  in  response  to  step 
torque,    6T  ,   or  step  command,  (A/3  +  B-^)     . 
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In  addition,  other  considerations  have  placed  the  following 
limits  on  K„  and  Kfi: 

K3  <    .40  (1.75) 

Kg   <    .021  (1.76) 

All  gains  (K, Kfi)  are  to  be  considered  positive. 

Thus  we  have  the  system,  and  the  performance  specifications 
to  be  met.     The  remainder  of  this  chapter  will  describe  briefly  the 
steps  to  be  followed  in  the  system  analysis. 

1.  4    Idealized  System 

Consider  a  system  in  which  the  servos  and  sensors  have  unity 
transfer  functions: 

Gl G5    =    1  (1.77) 

The  missile  with   "perfect"  sensors  and  servos  will  hereafter 
be  called  the  idealized  system.     The  idealized  system  will  be  con- 
sidered first,   and  then  the  effects  of  non-ideal  servos  and  sensors 
will  be  discussed.     It  is  sufficient  to  say  that  the  idealized  system  is  a 
reasonable  approximation  to  the  actual  system  if  it  can  be  shown  that 
the  servos  and  sensors  introduce  into  the  actual  system  additional 
modes  of  low  amplitude  which  are  far  removed  in  frequency  from  the 
natural  modes  of  the  idealized  system.     This  matter  will  be  discussed 
in  Chapter  6.     For  the  idealized  system  we  may  consider  the  aero1- 
dynamic  and  control  equations  in  the  following  form: 

(Ap  +B7  +M6)  =  p  -<l>  sin  a     +  $  cos  a  (1.78) 

C]3   +  E7  +N6  =  Tj?  (1.79) 

GS  +  F^'+Hj3  +L7  =  ?  (1.80) 

y  =  -g[  KX(A^+  B7+  M6)+  K2l  +  K3"]  (1.  81) 
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f*~[K4*+  K5J  +  Kg£]  (1.82) 

1.  5    Procedure  of  Investigation 

The  investigations  of  subsequent  chapters  will 

1.  Discuss  analytic  methods  of  selecting  gains  to  provide 
specified  performance. 

2.  Use  analogue  simulation  (guided  by  the  methods  discussed 
above)  to  select  a  suitable  set  of  gains. 

3.  Discuss  the  possibility  of  compensating  networks  within 
the  missile  to  eliminate  or  attenuate  the  effects  of  coupling  on  desired 
response. 
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CHAPTER  2 


ANALYTICAL  STUDY  OF  THE  COUPLED  SYSTEM 


2.1    General 

In  Section  1.  4  of  Chapter  1  are  shown  the  linearized  differential 
equations  of  the  idealized  system.     These  can  be  put  into  the  form  of  a 
"system  determinant"  by  gathering  all  of  the  variables  in  each  equation 
to  one  side,   the  other  side  being  zero.     If  this  be  done,   and  LaPlace 
transform  notation  be  adopted,  there  results  the  system  determinant: 


A  = 


0 

• 

• 

7 

6 

(S-A) 

cos  a 
o 

-sin  a 
o 

-B 

-M 

-C 

s 

0 

-E 

-N 

-H 

0 

(s-F) 

-L 

-G 

0 

0 

X 

2 

s 

0 

1 

-KXA 

-Y 

0 

(s-KjB) 

-KXM 

(2.1) 


where        X    =    (Kg  s      +    K.s  +  K4) 
Y    =    (K3s  +    K2) 


Expansion  of  this  determinant  yields  the  system  characteristic 
equation.     This  is  a  sixth  order  polynomial  in  s  in  whose  coefficients 
appear  a  considerably  involved  mixture  of  the  controller  gains. 
Appendix  B  gives  this  characteristic  equation  where  the  aerodynamic 
coefficients  for  a    =  24    have  been  entered. 

The  objectives  of  the  analytical  investigation  are: 

1.  To  determine  how  to  fix  upon  a  set  of  gains  K,  through 

14 


KR  which  produce  a  stable  system  and 

2.  To  determine  a  set  of  gains  which  will  provide  for  meeting 
the  specified  preformance  at  a  =  0  described  in  Chapter  1,  Section  1.  3 
and 

3.  To   select  these  gains  not  only  to  meet  objectives  1.    and  2.  , 
but  further  to  obtain  conditions  in  which  the  disturbance  in  roll  caused 
by  insertion  of  yaw  command  and  the  disturbance  in  yaw  casued  by  roll 
torque  disturbances  are  acceptably  low  and  of  reasonable  damping. 

Inspection  of  the  characteristic  equation  reveals  that  several 
inequalities  can  be  obtained  from  the  expressions  for  the  coefficients 
of  powers  of  s  which  must  be  satisfied  to  prevent  occurence  of  roots 
of  the  characteristic  equation  in  the  right  half  of  the  s  plane.     Given 
a  set  of  gains,   the  test  for  presence  of  sign  change  in  the  equation 
can  be  readily  performed. 

It  would  also  be  possible  to  apply  the  Routh  Criterion,   at  the 
expense  of  additional  labor,  for  a  given  set  of  gains,   to  obtain  further 
information  as  to  the  nature  of  the  roots. 

These,  and  other  procedures  that  might  be  derived  afford  little 
in  the  way  of  a  guide  to  selection  of  a  new  set  of  gains  to  test  when  one 
set  has  resulted  in  an  unstable  condition. 

There  is  the  possibility  of  forming  from  (2.1)  a  set  of  non- 
interacting  transfer  functions  from  which  a  block  diagram  of  the  system 
could  be  constructed.     With  manipulation  of  this  diagram  and  the  use  of 
root  locus  methods,   the  roots  for  a  given  set  of  gains  could  be  found. 
Unfortunately  there  would  be  at  least  three  inner  loops  of  such  a  diagram 
whose  roots  are  gain  dependent,   and  a  very  considerable  amount  of 
tedious  and  repetitious  graphical  work  would  be  involved,   even  if  each 
new  attempt  could  be  based  on  a  set  of  gains  choice  of  which  was 
guided  by  the  last  trial. 

It  is  apparent  that  it  would  be  desirable  to  have  a  more  orderly 
approach  in  which  less  futile  graphic  or  algebraic  work  is  expended  in 
arriving  at  the  conclusion  that  the  last  guess  was  unsatisfactory  before 
making  a  new,  more  or  less  intelligent,   guess. 
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The  object  of  the  following  section  is  to  discuss  a  procedure 
which,  while  not  avoiding  trial  and  error  and  not  avoiding  tedium,   at 
least  puts  to  good  use  the  experience  gained  by  past  results. 

This  procedure  begins  with  manipulations  of  the  characteristic 
equations  of  the  uncoupled  roll  and  yaw  subsystems.     These  equations 
are  of  low  order  and  quite  straight  forward  so  that  the  effect  of  gain 
changes  on  the  resultant  responses  are  readily  seen. 

2.  2    An  Algebraic  Approach  to  a  Stable  System  Characteristic 

Equation 

In  Chapter  3  is  derived  a  block  diagram  interpretation  of  the 
system  equations  which  results  in  a  new  form  of  the  system  determi1- 
nant  where  in  the  coupling  effects  and  the  subsystems  are  clearly 
depicted. 

When  (3.7)  oft  Chapter  3  is  reformed  ford£     =  0  we  have 


A    = 


Yaw  Subsystem 


1 

-F 

-F 

0 
0 


6 


-F 

1 
-F 

0 
0 


1 


-F 


12 


"F! 

1 

0 
0 


tv 


0 
0 
0 

-F 


10 


-F 


11 


(2.2) 


Roll  Subsystem 

Expansion  of  (2.  2)  will  reveal  a  result  which  is  the  product  of 
the  results  obtained  by  expanding  those  .sections  of  (2.  2)  labelled  Yaw 
Subsystem  and  Roll  Subsystem  as  if  they  were  themselves  determiL 
nants. 

When  the  block  diagram  Fig.    3.  7  of  Chapter  3  is  considered 
at   at     =  0  that  diagram  divides  into  two  new  diagrams  of  the  roll  and 
yaw  subsystems.     Using  the  methods  described  in  Chapter  3,   one  can 
obtain  the  subsystem  determinants  which  are  the  subdeterminants 
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within  (2.2). 

Then  if  we  let 
1 


10 


Fll       l 


s2  (s-F)  +  G(K6s2  +  K5  s  +  K4) 
s2  (s-F) 


(2.3) 


a,nd 


-F. 


-F 

1 
-F, 


13 


-F 


12 


-F, 


[  ^    +(-K1B-A-EK3)  v^2  +  (AKj  B-C-EK2  +  EAKg  -  BCK3  -  B^A)  ft 

i(s-K1B)(  ^Aj 


(CKLB  =  EAK2  +  EAKX  -  BCK2) 


§  (§-K1B)(s-A) 


(2.4) 


Employing  the  same  analog  with  Cramer's   Rule  as  is  used  in 

Chapter  3  we  can  express  the  outputs  due  (A/3  +  B7)      and   6T    in  this 

a     =  0    case  as 
o 


<Hs) 
6^TsT 


0  -F 

1  1 


10 


(2.5) 


(A/3  +  B(y(s) 
(A/3  +  BT(s) 


1 

-F        1 
1            0 

AFA 

"F3 

-F,       0 

BF 


A 


-F 


13 


1 


-F 
-F, 

1 


12 


(2.6) 


y  y 

When  (2.  5)  ,  (2.6)  are  expanded,   substitutions  made  for  the  defined 
transfer  functions,  and  simplification  made,   applying  a  step  input 
(A/3   +     B7)c(s)  =  ^       and    6T(s)  =  ^  . 
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*<S)=     -g ^—2 (2'7) 


(S    +  (GK6  -  F)  S"  +  K5  S  +  K4) 


-K    PA  [  Bs  -  Ej 

(A/3        +     B7   )(S)=    ^ rr— 

S  [  S    +  (-KjB  -  A  -  EK3)  S    + 


-K,  PB  [  S  (S-A)+  C] 

+ q ± n (2.8) 

S[  S    +  (-^B  -  A  -  EK3)S      + ] 


In  the  results  of  (2.  7)  and  (2.  8)  it  is  clear  that  the  modes  which 
exist  are  completely  specified  by  the  denominators  and  that  in  each 
case  these  are  3  in  number.     It  is  further  apparent  that  the  zeros  of 
both  closed  loop  functions  are  due  to  aerodynamic  coefficients. 

2   1 

Chu         has  derived  a  large  number  of  forms  of  "desired" 

characteristic  equations  on  the  assumption  that  the  performance  is 
dominated  by  a  pair  of  complex  roots  and  has  prepared  graphical  aids 
to  the  rapid  determination  of  the  pole-zero  configuration  of  the  closed 
loop  function  which  reflects  a  specified  time-domain  response  to  a 
step  input.     With  these  it  is  possible  to  rapidly  obtain  a  set  of  gains 
for  the  separate  subsystems  to  achieve  a  specified  response.     However, 
there  is  no  certainty  that  the  assumption  of  "2nd  order  dominance"  in 
the  subsystem  will  produce  a  set  of  gains  which  will  make  the  coupled 
system  stable. 

It  is  possible,    though,    to  assume  a  set  of  gains  and  rapidly 

determine  the  residue  coefficients  of  the  modes  that  exist  using  methods 

2    2 
described  by  Thaler.  By  this  method  a  check  can  be  quickly 

obtained  upon  the  ability  of  a  stated  set  of  gains  to  meet  the  specifi- 
cations at    a     =0.     A  map  can  then  be  prepared  of  the  locations  of  roots 
of  the  characteristic  equations  of  the  subsystems  which  do  satisfy  the 
a  =  0. specifications. 

Having  obtained  a  starting  set  of  values  of  the  gains,   these  can 
be  entered  into  the  coupled  system  characteristic  equation  (Appendix  B) 
and  that  equation  factored,   the  roots  plotted  with  the  roots  of  the 
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characteristic  equations  of  the  uncoupled  subsystems,   and  the  trend 
noted. 

A  second  trial  for  the  roots  can  be  guided  by  the  concept  that 
increased  stability  of  the  coupled  system  will  be  accompanied  by 
increased  "sluggishness"  in  the    a     =0  modes  of  the  subsystems.     In 
the  actual  carrying  out  of  this  process  it  is  also  noted  that  transient 
oscillating  frequencies  of  a     =  0  modes  tend  higher  from  the  first 
choice  of  subsystem  root  locations.     By  following  a  path  on  which  the 
complex  pair  of  roots  of  each  subsystem  increase  in  natural  frequency 
but  remain  at  constant  damping  ratio  the  real  root  moves  toward  the 
origin,   and  convergence  to  stability  of  the  coupled  system  characteristic 
equation  is  apparently  more  rapid  than  on  other  paths. 

Thus  there  does  exist  one  orderly,  albeit  tedious,  process  for 
converging  upon  a  set  of  gains  which  produces  stability  of  the  coupled 
system  and  yet  allows  the  subsystems  to  meet  the    a    =  0  specifications. 

"Bold  spotting"  allowed  convergence  to  a  stable  set  of  gains  in 
three  tries  using  these  concepts. 

Somewhat  less  tedious  labor  is  involved  if  one  does  not  observe 
the  constraint  of  specifications  too  rigorously,   merely  working  with  the 
subsystem  characteristic  equations  to  produce  stability  in  the  coupled 
system  characteristic  equation  without  too  much  regard  for  specifications, 
then  trimming  to  achieve  the  specifications  at  a-    =0. 

Such  a  procedure  is,   then,    capable  of  satisfying  objective  (1)  of 
Section  2.1.     It  is  quite  apparent  that  repetitive  factoring  of  the  coupled 
system;  characteristic  equation  as  would  be  necessary  with  this  (and 
most  other)  methods  to  reach  the  other  objectives  is  simply  beyond 
reason  without  automatic  computing  facilities. 

The  conclusion  is  therefore  drawn  that  this  method  does  set  forth 
a  concept  which  is  useful  in  guiding  one  to  a  stable  set  of  gains  for  an 
analogue  simulation  of  the  problem,  but  is  otherwise  of  academic  interest, 
and  totally  devoid  of  practicality  in  achieving  the  other  objectives  of 
Section  2. 1. 

For  this  reason  all  work  in  "trimming  up"  the  system  was  done  by 
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analogue  simulation  as  discussed  in  Chapter  4. 

2.  3    Steady  State  Analysis  of  the  System 

In  the  preceding  section  it  was  concluded  that  setting  the  gains 
in  the  system  to  meet  all  the  stated  objectives  is  not  a  practical  task  for 
analytical  study,    that  rather  it  is  an  undertaking  more  appropriate  to 
the  analogue  computer.     Considerable  useful  information  for  a  com- 
puter study  can  be  gained  by  analyzing  the  features  of  the  variables 
considered  as  outputs  in  the  steady  state. 

The  methods  described  in  Chapter  3  for  solving  the  equations 

will  be  used  here.     The  uncoupled  mode  at  a    =  0  will  be  dealt  with 

^  o 

first. 


1.     Uncoupled  Roll  Subsystem 

In  Section  2.  2  it  was 
toristic  equation  is  given  by 


In  Section  2.  2  it  was  shown  that  at    a    =0  the  roll  system  charac- 

o  J 
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-F 


-F 


10 


11 
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S2(s-F)  +  G  (Kg  £2  +  K5s  +  K4) 
a2  (,s  -  F) 


(2.3) 


It  was  also  shown  that 


0      -F 


1         1 


4-  <s>- 

6T  A 


10 


10 


(2.5) 


Since 


10    , 


applicable.     For       6T(s)  = 


is  analytic  in  the  RHP,   the  final  value  theorem  is 
K 


R 


Lim  </>(t)  =  Lim  [  S<KS)J 
t-  cO  s  -  0 


a (s  -  f) 


S  +  (GK, 


F)s    +  GK5s  +GK4 


32  (,s  -  F) 


=  0 


(2.9) 
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Thus  the  roll  system  ultimately  recovers  from  a  roll  step 
disturbance. 

2.     Uncoupled  Yaw  Subsystem 

From  (2.  4)  the  characteristic  equation  of  the  uncoupled  yaw 
subsystem  is 

[  s3  +(-1^  -  A  -  EK3)  s2  +(AK1B  -  EK2  +  EAK3  -  BCK3  -  BK}A)s 
Ay  =  s  (s   -  KL  B)(s  -  A) 

(C^  B  -  EAK2  +  EAKX   -  BCK2)] 
+  s  (s  -  KJBRi  -  A)  (2'4) 

In  (2.8)  for  a  step  input  (A /3+  B7)    (s)  =  — 

-  PAK    [  Bs  -  E] 
(A/3+BT)(s)    =    3- 


s  [  s     + (BC  -  AE)(K,  +  K2); 


-PBK    [  s2  -  As  +  CJ 

+    . i (2.5) 

s  [    s     + (BC-  AE)(K1  +  K2)J 


Then  since  for  a  stable  choice  of  gains  this  function  is 
analytic  in  the  RHP  the  Final  Value  Theorgm  applies  and 

Lim  (Aj3+  B7)(t)  =  Lim  [  s  (Aj3+  Bt)(s)] 
t-OO  s-*0 

-PK,  (BC  -  AE)  K, 

P( 1 ) 


K1  +  K2(BC-AE)  K1+K2  (2   n) 


Applying  the  Cramer's  Rule  method  of  solution  to  the  yaw 
subdeterminant  in  (2.  2) 

•  -Fa   (-F6  -F3F5)      [  Es(s-A)+BC]     s(s_A)(s_KiB) 


-(s) 


(A/3  +  B7)c—  Ay  Ay  (212) 
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or 


± 


(A/3+B7) 


(s)  = 


-K    [  Es  +  (BC  -  AE)J 


[s3  + 


(BC  -  AEXKj^  +  K2)] 


(2.12a) 


Using  (Aj3  +  B7)   (s)  =  —  again,   and  applying  the  Final  Value 


Theorem 


"Kl  p 

Lim    ty  (t)  =   Lim  s    [  $  (s)]  =  „    ,  g 

t-00  s  -*  0  Kl+K2 


(2.13) 


Thus  it  is  apparent  that  both  (A/3  +  BY)(t)  and  i|i(t)  come  to  a 
constant  non  zero  value  in  steady  state  and  that  4<(t)  increases  without 
limit   .     Their  senses  and  magnitudes  are  easily  calculated  for  com- 
parison with  computer  results. 


3.     Coupled  System 

In  performing  the  steady  state  analysis  at  a    >  0  it  is  convenient 
to  let  s  -*■  0  in  the  system  determinant  (3.  7)  at  the  outset  since  that 
determinant  appears  in  the  denominator  of  the  Cramer'  s  Rule  solution 
for  each  variable  of  interest.     In  the  following  statements,   the  fact  that 
this  has  been  done  will  be  signified  by   [  A]         where 


Lim  A 
s  -  0 


K, 
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-E     0 


i-  ° 
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F 
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F 
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F 

-3tA'ss 


(2.14) 


Other  notation  is  as  explained  in  Chapter  3.     It  is  assumed 
throughout  that  a  stable  situation  exists  so  that  the  Final  Value  Theorem 
is  applicable. 
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a.     (/>  due  to  roll  disturbances  &T(s)  =  


</>(s)  =  S     Arrs  (2.15) 


S^5-[  A  o       2    r    a 

T2  L      rrJ  ss  R  s     [  A  ^  J 


Lim  </>(t)  =Lim  s  [  <j>(s)]  =    -^ =    rr    ss  (2.15a) 

t— o  b-0  -^-[A]gs  lA^ss 


Evaluation  of  the  determinants  shows  </>(oa)  =  0  when  a  roll 
disturbance  is  applied. 

p 

b.  §  due  to  yaw  command  (A/3  +B7)   (s)  =  — 

s     yr  s 
<Ms)=  — (2.16) 

s^r[  A      ] 

I  L      yr    ss 

Lim0(t)=Lim     s-[^(s)]  =  — ~ (2.16a) 

t^co  s-*0  -3[A]SS 

s 

Again  evaluation  of  the  determinants  shows  <£(«<>)  =  0  when  a  yaw 
command  is  applied. 

•  P 

c.  ip  due  to  a  yaw  command  (Aj3+ B7)   (s)  =  — 

c  s 

A  P 


i|;  (s)=    W    S  (2.17) 


•   .  .  .    • 


P  s3  I  AvvJ 
Lim  <\>  (t)  =    Lim  s  [  cJj(s)]  =  ^~ 

t-oo  s^O  t  A^  ss 


ss 


(2.17a) 


Evaluation  of  the  determinants  shows  that  \\>  goes  to  the  constant 
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-K    P 

[+(01        =  I? — T7-  (2.17b) 

LT       Jss       Kn  cos  a     +  K0 
1  o  I 


and  that  i|j  increases  without  limit  as  t-» 


d.     ^  due  to  roll  disturbance    6T(s)  =  — 


^A 
Us)=- &-  (2.18) 

A 


.               Rs3[A       ] 
Lim  4j  (t)=  Lim  s  [  ^  (s)]= r$     ss  (2.18a) 

t  -*o©  s  -*  0  "■      '  ss 


Evaluation  of  the  determinant  shows  that  vjj  (°°)  =  0. 

Further  evaluation  of  the  integral  shows  that  \\t  (p°)  =  0. 

This  information  on  steady  state  values  is  of  use  in  correlating 
the  analogue  computer  results  in  Chapter  4  with  the  physical  operation 
of  the  system. 
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CHAPTER  3 
COMPENSATION 


3.1    General  Approach 

The  general  problem  of  preventing  interaction  between  con- 
trolled quantities  in  multiple-loop  systems  has  been  discussed  by 

3   13    2  2   1 

Tsien,    *     Truxal   '     and  others.     Among  these  Chu   *     has  developed 

a  generalized  theory  of  multiple-loop  systems  including  techniques  for 

the  determination  of  compensation  necessary  to  prevent  interaction. 

Chu  employs  methods  of  manipulation  with  the  system  determinant 

with  the  aid  of  block  diagrams.     The  methods  of  Tsien  and  Chu  are 

quite  parallel  in  theory,   the  choice  between  them  being  a  subjective 

matter.     Chu'  s  method  has  been  chosen  for  application  to  this 

problem  since  that  technique  is  more  familiar  to  the  authors. 

A  brief  description  of  the  basic  elements  of  Chu.'  s  approach 
to  mulitple-loop  control  systems  is  included  in  Appendix  C  and  no 
further  detailed  explanation  will  be  given  here. 

The  problem  of  analyzing  the  effects  of  interactions  in  this 
system  is  essentially  one  of  maintaining  contact  with  the  physical 
hardware  involved.  Particularly,  a  block  diagram  possesses  the 
general  appeal  of  being  a  reasonably  tangible  representation  of  the 
relations  among  the  physical  variables  of  a  system,  provided  that 
diagram  does  not  endeavor  to  lump  too  many  quantities  and  inter- 
relations into  a  single  complicated  block  transfer  function. 

A  diagram  which,  it  is  believed,  will  provide  a  reasonably 
physical  view  of  the  problem  will  first  be  derived.  It  employs  the 
most  basic  set  of  transfer  functions  obtainable  from  the  linearized 
equations. 
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Consider  the  set  of  linearized  equations  in  the  form  of  a 
determinant: 
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and 
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(K6  s     +  K5s  +  K4) 
(K2  +  K3s) 


From  each  equation  represented  by  a  row  of  the  determinant 
a  single  variable  is  selected  to  be  described  as  a  function  of  all  the 
others.     For  manipulative  ease,   a  transfer  function  describing  each 
relation  is  defined: 
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Equation  (I) 
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(3.2a) 
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i 


sin  a 


(s)   = 


OT    '=     F2<S> 


(3.2b) 


^    =    OT     ='    F3<S> 


(3.2c) 


T(s)  =  r^A)  s  F4(s) 

0 


(3.  2d) 


The  theorem  of  superposition  suggests  representation  of  |3(s) 
by  the  elemental  block  diagram. 
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Fig.    3.1 
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Equation  (II) 
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(3.4) 


hB)  -i^)  -=  Fs(s) 


(3.4a) 


T 


(s)  = 


T^fT 


=    FAs) 


(3.4b) 


G 


4<s>  =  T^f)       Fio(s 


(3.4c) 


28 


M*L 


iLsl 


M§L 


10 


<?(s) 


Fig.   3. 3 
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Equation  (IV) 
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It  will  be  noted  that  provision  has  been  made  for  the  insertion 
of  a  yaw  command  signal  which  subtracts  from  the  sensed  side 
acceleration  (Aj3  +  B7  +  M6)  as  indicated  by  Figure  1.  3  of  Chapter  1. 
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Fig.    3.5 


The  interrelations  existing  among  the  system  variables  may 
then  be  described  by  a  block  diagram  as  shown  in  Figure  3.  6. 
Provision  of  a  node  for  introduction  of  a  roll  disturbance  torque  6^ 
will  be  noted. 
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Fig.    3.6    System  Block  Diagram 
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The  method  of  Chu  requires  translation  of  the  block  diagram  of 
a  system  such  as  depicted  in  Figure  3.6  into  a   "Standardized  Block 
Diagram".     In  this  "standardized"  diagram  certain  manipulations  to 
Figure  3.  6  are  made  to  satisfy  a  set  of  conventions  which  allow  a 
symmetry  to  exist  between  the  "standardized"  diagram  and  the  system 
determinant  and  which  allow  the  writing  of  the  system  determinant 
by  inspection.     These  are  described  in  Appendix  C.     When  the 
necessary  adjustments  are  made  the  "standardized"  block  diagram 
takes  on  the  appearance  of  Figure  3.7. 

In  Figure  3.  7  it  will  be  observed  that  the  manner  of  redrawing 
Figure  3.  6  to  produce  the  "standardized  block  diagram"  has  resulted 
in  the  division  of  the  coupled  system  into  regions  which  can  be  referred 
to  as  the  roll  and  yaw  subsystems.     It  is  readily  apparent  that  Fq  and 
Fq  represent  the  coupling  effects  of  yaw  into  roll  while  F?,   F.,  F7  and 
F,.  clearly  show  the  coupling  effects  of  roll  into  yaw.     Inspection  of 
the  transfer  functions  represented  by  F?,  F.,  F?,  Ffi,  Fq>  F,,  immedi- 
ately show  that  at  a    =  0  the  system  separates  into  two  subsystems. 

By  application  of  the  rules  (Appendix  C)  for  writing  the  system 
determinant  from  the  "standardized  block  diagram",  that  determinant 
becomes: 

Yaw  Subsystem 


1 


-F, 


-F. 


-F, 


-F 

1 

-P, 


13 


-F 


12 


Roll  to  Yaw  Couplings 
-F, 


-F, 


0 

0 

-P. 


14 


-FP 


-F 


4 


0 


-F, 


-F 


10 


-F 


11 


(3.7) 


Yaw  to  Roll  Couplings 


Roll  Subsystem 
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Equating  A   to  zero  and  evaluating  one  obtains  the  system 
characteristic  equation.     Immediate  proof  that  above  represents  the 
original  determinant  of  the  system  can  be  obtained  by  reading  each 
row  of  A  with  the  appropriate  substitutions. 

As  a  consequence  of  the  way  selected  for  redrawing  Figure 
3.6  to  produce  Figure  3.7,  A   divides  into  regions  illustrative  of  the 

subsystems  and  the  intercouplings  as  noted  above.     The  close 

3   1 
parallel  between  this  determinant  and  Tsien'  s         matrix  are 

apparent.     Also  the  clarity  with  which  interactions  are  displayed  in 

Figure  3.  7  and  in  A  would  seem  to  some  degree  to  refute  one  of  the 

3    3 
objections  to  the  block  diagram  approach   "    . 

It  is  true  that  considerable  effort  (although  simple  and  direct) 
has  been  put  into  construction  of  this  form  of  A  ,  but  this  has  placed 
the  designer  in  the  position  of  having  a  graphic  representation  in 
Figure  3.  7  of  Cramer'  s  Rule.     Since  each  node  output  of  Figure  3.  7 
represents  one  of  the  variables  of  interest  in  the  problem  it  is  merely 
necessary  to  follow  Cramer'  s  Rule  to  assess  the  output  at  any  node 
due  to  input  at  any  node. 

3.  2    Objectives  and  Procedures 

The  objectives  of  complete  compensation  will  be  assumed  to 
be: 

1.  isolation  of  </>     from  (A/3  +  By) 

2.  isolation  of  \\)  from    6™ 

The  introduction  of  any  compensation  will  change  A  ,   and  in 
this  will  clearly  affect  all  outputs.     This  change  should  not  produce 
serious  effect  in  the  desired  responses,     (i.  e.     Should  not  produce 
an  unacceptable  array  of  closed  loop  poles). 

From  the  foregoing  statements  regarding  the  association 
between  Cramer'  s  Rule  and  Figure  3.7  and  A   it  is  clear  that  it  is 
possible  to  solve  for   <f>  (s)  as  a  function  of  (A/3  +  By)     as  follows: 
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W*+BT)    M-^X.  (3.8) 


c 


Where  A        has  been  formed  from  A  by  replacing  the  IV  (or  </>) 

column  of  A  by  a  column  having  zeros  in  all  row  elements  except  the 

first  at  which  unity  appears  (the  entrance  function  F     and  the  yaw 

a  • 

command  signal  transform  having  been  factored  out  of  the  <f>  column  of 


A      ). 


Likewise 


The  desired  response  in  yaw  is  given  by 

jj-i-^s,  ,  _^0L  (3.10) 

And  the  normal  (existing,   though  not  particularly  desired) 
reaction  of  <£   to  6™  is 


Provided  the  missile  can  be  prevented  from  yawing  due  to  roll 
disturbance,  no  side  acceleration  will  occur. 

Then  the  desired  goals  of  compensation  are  achieved  if 

1.  A        =   A         =  0 

yr  ry 


TTn  — 77T  (s)    is  not  altered  by  accomplishment 
(Aj3  +  B7)c 

of  (1)  except  in  the  sense  that  A  is  changed,  and 

that  this  change  is  not  deleterious. 
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It  is  of  interest  to  inspect  both  A         and    A 

*-  yr  ry 


yr 


-F, 


-F. 


-F, 


0 


— ¥*-   —  -F-~ 
^13       ~2 


1  -F5  f 

-Fl  l  f 

0  -F8  f 

0  0  |) 


14 


-F, 


-F 


4 


-F 


10 


(3.13) 
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-F,         0 
o 


-F  1 

*8         l 


-Q F 


-F 


14 


-FP 


-Fn      -F 


-F 


10 


11 


(3.14) 


As  the  dashed  lines  indicate,   the  expansion  of  A         and  A 

f  ry  yr 

by  cofactors  pivots  on  the  circled  element  and  suppresses  the 
indicated  rows  and  columns.     Thus 

1.  a  compensator  placed  at  rl,   c  IV  will  not  appear  in 
A 

yr 

2.  a  compensator  placed  at  r  V,   ell  will  not  appear  in 
A 

ry 

Also  of  interest  are  A        and  A       : 

rr  yy 
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(3.15) 


yy 


-0- 


-F, 
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-F, 


12 


-Ff 


14 


0  -Fr 


1 

"F2 

"F4 

-F 
*8 

1 

"F10 

0 

"Fll 

1 

(3.16) 


From  these  it  is  noted  that  a  compensator  placed  at  rl,   cIV  does 

not  appear  in  A        and  a  compensator  placed  at  rV,   ell  does  hot  appear 
in  A 

yy 

Thus  one  concludes  that  compensation  at  rl,   cIV  can  possibly 
be  used  to  make  A        =0  and  a  compensator  at  rV,   ell  can  possibly  be 
used  to  make  A        =0  and  that  neither  compensator  will  appear  in 

A        and  A 

yy  rr 

From  the  way  in  which  the  system  determinant  is  constructed 
from  the  "standardized  block  diagram"  it  is  apparent  that  the  rl,   cIV 
element  of  A   constitutes  a  feedback  from  node  IV  to  node  I  in  Figure 
3.7.     Inspection  of  the  diagram  reveals  that  use  of  this  path  implies 
that  <j)   is  available  and  that  modified  </>   can  be  added  algebraically  to 
7  .     Likewise  rV,   ell  of  A     constitutes  feed  forward  from  node  II  to 
node  V.      ijj   i-s  physically  tangible.     The  same  question  of  ability  to 
mix  the  modified   41    signal  with   6    arises  here.     If  the  compensating 
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networks  required  are  realizable  either  exactly  or  approximately  this 
question  of  mixing  can  be  studied  further. 

Installation  of  compensation  functions  -  V(s)  at  rV,   ell  and  W(s) 

at  rl,   cIV  in  A  ,  where,   conveniently,  no  element  previously  existed, 

results  in  the  solutions,   when  A        and  A        are  equated  to  zero: 

yr  ry 

-Ls     +  (LA  -  BH)s  +  (HE  -  LC)  cos  aQ 
V(s)  =  (GE-LN)s  +  (-AGE-LMC-HNB+LAN  +  HME+GCB)  <3-17) 


,r/  1  612s2  +  143s  -  146500  /Q   1Qv 

*a   =24°  =   ^Os~T42200) (3*18) 

o 

.  1  _    133.5s2  +  47.8s  -  1548  .„   1Q, 

KSJ     =12°     -[  (100>300)s  +  37,570]  ^.ih; 

sino(D(s-A)[(CG-HMK1)s+K1(-BCG-LAN-HME+LCM+AGE] 


W(s)  = 


(s-K1B)[(GE-LN)s+(-AGE-LMC-NHB+LAN+HME+GCB)] 

s  (s-F)[  Ns+(CM+K,ME-NA-BNK1)] 

+    

(s-K1B)[(GE-LN)s+(-AGE-LMC-NHB+LAN+HME+GCB)] 

(3.20) 

Unfortunately  neither  of  these  networks  is  realizable  with 

3.4 
passive  circuit  elements  exclusively. 

Upon  examination  of  V(s)  it  is  found  that  it  can  be  divided  into 

2 


v(s; 


1  _  6.71s  1.57s  1600  ,„  . 

5J      _Q.o=    "    s+463       "  s+463     +  "sT463  l^.iHa) 


kv 

6.71s             1.57s          1600 
s+463           s+463       '   s+46  3 

1 

1.  335xl0"3s2       .478xl0"3s      15.5xl0-3 

-19° 

o 

s  +  .376               s  +  .  376        '     s  +  .  376 

„.   1  _      1.  335x10      s         .478x10      s      15.  5x10  ,„   1Q   . 

(SJ      _19o  =  "      s  +  ,376  s  +  .376       +    s  +  .  376  U.  iaa) 
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An  approximation  using  the  last  two  terms  of    V(s)  is 

J24° 
employed  in  the  analog  computer  simulation  with  the  result  that 

while  addition  of  this  approximation  to  V(s)  does  not  create  an 

unstable  system,  neither  does  it  attenuate  the  roll  disturbance  due 

to  a  yaw  command  input.     Thus  the  indication  is  that  the  first  term 

is  of  dominant  importance. 

Thus  hope  of  a  simple  approximation  to  V(s)  which  would 

attenuate  the  unwanted  roll  response  is  lost.     By  arrangement  of 

the  analog  simulator  to  provide  an  approximation  to  the  first  term  of 

V(s)  above,   about  30%  decrease  in  the  extremes  of  roll  response 

24 
to  yaw  command  is  achieved,  indicating  that  at  least  partial  decoupling 

can  be  achieved.     Addition  of  the  second  term  as  above  does  not 

show  any  appreciable  improvement  over  the  first  term  alone.     At 

this  point  available  facilities  on  the  computer  are  exhausted.     The 

actual  computer  techniques  employed  are  discussed  in  Chapter  4. 

Figure  5.  7  shows  the  improvement  achieved. 

Since,  if  isolation  is  to  be  achieved,   simulation  of  the  first 

term  of  V(sM         is  of  paramount  importance,  it  would  appear  that  a 

J24°  # 

tandem  pair  of  derivatives  on  the  available  \\t  sensor  output  is 

required.     In  this  connection  it  is  worthy  of  note  that  the  system  can 

be  driven  unstable  quite  easily  when  more  gain  than  specified  is  used 

in  the  double  derivative  term  of  V(s). 

The  problem  of  introducing  the  modification  to  mix  with  6  has 

3   5 
not  been  examined  here.     Available  literature   *     does  not  indicate 

any  hope  of  introducing  the  compensating  signal  other  than  to 

differentiate  it  and  mix  it  with  6     rather  than  6.     That  this  procedure 

works  has  been  verified  on  the  analog  computer. 

It  will  be  noted  in  Figures  5.  8  and  5.  9  that  the  desired 

« 
responses  to  a  yaw  command,  (A/3  +    B7  +  M6)  and  i|j  ,  have  not  been 

adversely  affected  by  the  approximate  compensation  employed. 

Relative  to  the  possibility  of  using  an  expanded  version  of 
W(s)  similar  to  that  employed  to  approximate  V(s)  ,  it  will  be  noted 
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that  W(s)  is  dependent  on  K,,   and  thus  the  process  of  adjusting  gains 
would  involve  adjustment  of  the  compensator.     This  would  result  in  a 
"tail-chase"  operation  and  is  not  felt  to  be  practicable.     Computer 
facilities  available  do  not  permit  examination  of  the  effects  of 
approximate  simulation  of  W(s). 

Since  the  desired  approximate  compensation  for  a    =12     and 

a    =  24     are  widely  different,   it  is  apparent  that  if  the  compensation 

for  a    =24     is  installed  permanently,   the  compensation  for  a    =12 

will  not  be  proper,   even  in  an  approximate  sense.     Also  any 

compensation  left  in  permanently  will  have  the  effect  of  creating 

coupling  at  a     =0.     Then  some  form  of  logical  operation  is  necessary. 

It  is  suggested  that,   since  the  a    =0  and  a    =12     responses  are  not 

greatly  different  (See  Figures  5.1  and  5.  2),   the  compensator  for 

a    =  24     alone  be  installed  and  that  this  be  switched  in  and  out. 
o 

Inasmuch  as  maneuvers  which  would  be  likely  to  develop 

extreme  angles  of  attack  would  probably  be  presaged  by  some  violent 

excursion  of  the  seeker  tracking  line  relative  to  the  airframe 

longitudinal  axis,   it  seems  conceivable  that  this  information  could  be 

used  as  a  basis  for  the  necessary  compensation  switching.     This 

possibility  has  not  been  examined  in  detail. 

In  conclusion,  a  method  has  been  described  for  discussing 
compensation  in  a  clear  manner  which  involves  explicit  relation  among 
the  measured  quantities  available.     A  compensation  which  will 
theoretically  completely  isolate  roll  output  from  yaw  command  has 
been  found  and  has  been  approximately  verified  on  the  analog 
simulator. 
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CHAPTER  4 


COMPUTER  SIMULATION 


4. 1    Description  of  the  Computer  Setup 

The  missile  system  is  simulated  on  a  Reeves  Electronic 
Analog  Computer  (REAC).     Outputs  are  recorded  on  a  Sanborn  4 
channel  recorder.     The  idealized  system  is  instrumented,  and  can  be 
represented  by  equations  (4.1)  through  (4.5)  below.     These  are 
rearrangements  of  equation  (1.78)  through  (1.  82),  and  provide  5 
equations  in  5  unknowns  with  each  equation  written  to  specify  the 
highest  derivative  present  of  one  of  the  variables.     The  functional 
layout  of  the  equations  is  shown  in  Figure  4. 1.  * 

*l     =     G(6+6T)   +F(j)   +    Hj3   +  LT  (4.1) 

6     =     -(K4</>  +  K5J  +  K6V)  (4.2) 


•  •  • 

$     =     (Aj3+  B7+M6)   +  sin  a    </>    -  cos  a     \\i  (4.3) 

tj?    -     Cj3  +  E  7  +  N6  (4.  4) 

7     =     Kx((A/3+ B7 +M6)-(A(3  +  B7 +M6)c] 

+    K04>  +  K  *if  (4.  5) 


cFigure  D.  1,  Appendix  D  shows  the  actual  scaled  layout  used. 
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The  system  of  Figure  4.1  is  the  one  on  which  all  gain  variation 
and  adjustment  tests  are  run.     In  addition,   the  compensator  network  of 
Figure  4.  2,   representing  the  approximate  compensator, 


f  (s)  =V(s)] 


* 


24^ 


6.71s 
s+463 


1.57s 
s+463 


(4.6) 


as  developed  in  Section  3.2  is  used  in  testing  the  effectiveness  of 
decoupling  compensation  at    a    >   0. 

Some  modifications  of  the  layout  of  Appendix  D  are  necessary 
to  permit  the  insertion  of  the  scaled  compensator  as  a  transfer  function 
from  ^    to  6.     The  modifications  are  minor,  however,  and  the  system 

is  not  shown  in  its  modified  form.     The  approximate  derivative  network 

4  1 
is  from  Johnson,  Analog  Computer  Techniques.    * 


Note:    RC  = 
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Fig.   4.  2     Analogue  Computer  Simulation  of  Decoupling 

Compensator 
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4.  2    Gain  Adjustment  Techniques 

The  following  gains  found  by  the  methods  of  Section  2.  2    to 
stabilize  the  system  are  used  as  a  starting  set  in  adjustment  to  meet 
specifications. 

K       =        43.8  (4.7) 

K2  =  6.18  (4.8) 

K3  =  .357  (4.9) 

K4  =  27.9  (4.10) 

K5  =  2.33  (4.11) 

Kg  =  .01985  (4.12) 

Roll  adjustment  is  done  by  recording  and  observing  </>;  yaw 
adjustments  made  on  the  recorded  step  response  of  (A0+  By  +  M6). 

By  observation  of  the  recorded  responses,   cyclical  gain 
adjustments  in  roll  and  yaw,  and  at  the  3  pitch  angles  selected  are 
made.     It  is  found  that  K„  and  Kfi  are  to  be  set  at  or  near  their 
maxima.     Next  adjustments  of  K.  and  K,-  in  roll,  K?  and  K,  in  yaw, 
are  made  xto  give  specified  stable  operation  of  the  subsystems  with 
suitable  damping  of  the  coupled  modes.     Subsystem  frequencies  (and 
thus  rise  time)  are  dependent  on  K,  and  K   •  damping  p,n  Kp  and  K_  . 
These  are  the  rough  guides  used  in  obtaining  the  final  gains  listed  in 
Section  5.1.     Section  5.2  sjhows  system  response  to  step  6T  and  to 
step  (A/3  +  By+  M6)     at  the  gains  selected.  * 

4.  3    Compensator  Simulation 

A  scaled  setup**  of  the  compensator  of  Figure  4.  2  is  used  in 
testing  the  possibility  of  eliminating  or  attenuating  the  effects  of  a  yaw 
command  on  roll  output. 


^Figure  D.  3,   Appendix  D  shows  correlation  between  analytically 
derived  response  of  the  uncoupled  mode  and  computer  response  for 
a  common  set  of  gains. 

**See  Figure  D.  2,  Appendix  D. 
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Responses  of  the  system  with  and  without  compensator  are 
compared  in  Figure  5.7  through  5.9. 
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CHAPTER  5 


RESULTS 


5. 1    Gains  Selected 


Final  gains  selected  are: 

Kx      =        32.0  (5.1) 

K2     =  6.40  (5.2) 

K3     =  .357  (5.3) 

K4     =        12.85  (5.4) 

K5     =  2.33  (5.5) 


K6     =  .0199  (5.6) 


5.  2    System  Step  Responses 

Figures  5.1  through  5.9  are  responses  of  <f>,  (A/3+By  +  M6), 
and  i)j  to  step  inputs: 


6T     =        1  degree  (5. 7) 


(A/3+B7+M6)       =    1  deg/sec  (5.8) 
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5.  3    System  Response  Parameters 

Reduction  of  response  data  from  Figures  5.1  through  5.6  gives 
the  uncompensated  system  response  parameters  show  in  Table  5.1. 


Roll 

Yaw 

deg 

roll 
sec. 

rmax 
deg. 

Tyaw 
sec 

(A/3  +B7  +  M6) 

max 

deg 
sec 

Input 

0 

.13 

.  56 

0 

6T 

=  1° 

12 

.14 

.61 

.06 

24 

.  34 

.61 

.25 

0 

y 

0 

.44 

.84 

(A/3  +B7)C 

_  1  deg 
sec 

12 

.48 

.  35 

.87 

24 

.57 

.23 

.94 

Table  5.  1 


Response  Parameters  of  the  Uncompensated  System 

Table  5.  2  is  reduced  data  from  Figures  5.7  through  5.  9  and 
gives  parameters  of  the  compensated  system  at     a     =  24°, 


Troll 

ymax 
deg 

vaw 
sec 

(A/3+B7  +M6) 

max 

sec 

deg/sec 

.34 

.43 

.24 

.91 

Table  5.  2 


Response  Parameters  of  the  System  as 
Compensated  at  a    -  24 


*Roll  and  yaw  subsystem  time  constants  as  defined  in  Section  1.  3. 
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CHAPTER  6 


CONCLUSIONS 


6.1    Achievement  of  Specifications 

The  uncoupled  response  specifications  of  Section  1.  3  can  be 
readily  achieved.     At  or     >  0,   coupling  exists  which  must  be  elimin- 
ated or  damped  to  maintain  stable,   specified  performance.     Stable 
coupled  response  is  obtained  at  the  expense  of  more  sluggish 
performance  in  the  uncoupled  mode. 

Satisfactory  response  in  both  coupled  and  uncoupled  modes 
may  be  obtained,   but  requires  computer  simulation  for  practical 
investigation  and  setting  of  system;parameters. 

6.  2    Elimination  of  Coupled  Response  by  Compensation 

Decoupling  by  compensation  is  theoretically  possible.     However, 
the  compensating  devices  cannot  be  achieved  with  purely  passive 
networks.     Also,    some  physical  changes  within  the  control  systems 
are  necessary  to  provide  points  at  which  compensating  signals  may  be 
introduced  and  extracted.     Third,   since  compensation  is  a  function  of 
angle  of  attack,   there  should  be  some  manner  of  varying  the  compensator 
characteristics  with  angle  of  attack. 

Since  some  success  has  been  achieved  with  approximations  of 
theoretical  compensators,   a  more  intensive  investigation  of  decoupling 
compensators  might  evolve  a  compensator  which  is  both  effective  and 
physically  practical. 

6.  3    Effect  of  Non-Ideal  Servos  and  Sensors 

Since  the  bandwidths  of  sensing  devices  and  control  surface 
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servos  commonly  employed  in  missiles  of  this  type  are  of  the  order  of 
40  to  50  cps  it  is  clear  that  the  factors  of  (B.  25)  are  quite  considerably 
removed  from  the  poles  (s-A)  and  (s-F)  and  from  the  origin. 

The  results  of  factoring  the  ideal  system  characteristic  equation 
indicate  that  its  roots  are  not  fG^r  removed  from  the  poles  thereof  for  the 
gains  selected  by  computer  study  in  Chapter  4.     Then  by  inference,    such 
additional  roots  as  are  introduced  in  the  non-ideal  system  by  the 
presence  of  control  surface  servos  and  sensing  devices  would  be  located 
relatively  close  to  the  poles  introduced  by  those  devices,   and  the 
residues  associated  with  these  roots  would  be  small. 

Therefore,   mathematical  representation  by  the  idealized  system 
is  considered  correct  and  no  gross  performance  characteristics  are 
thereby  over  looked.     Additionally,   as  noted  elsewhere,  available 
analogue  computer  facilities  do  not  permit  simulation  of  the  servos  and 
sensing  devices.     Conversation  with  Professor  W.  L.    Markey  of  M.  I.  T. 
relative  to  the  relation  of  additional  modes  introduced  by  non-ideal 
instruments  to  the  basic  modes  of  such  systems,  (especially  when  these 
instruments  have  wide  bandwidths  compared  to  the  natural  modes  of 
the  system)  would  seem  to  comfirm  that  this  simplifying  action  is  in 
keeping  with  common  engineering  practice. 
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APPENDIX  A 


AERODYNAMIC  COEFFICIENTS 


Values  of  the  aerodynamic  coefficients  defined  by  equations 
(1.  53)  through  (1.62)  and  evaluated  at  the  selected  values  of  a 
considered  in  this  thesis  are  shown  below  in  Table  A.  1.     The  values 
were  obtained  by  Convair  (Pomona)  engineering  personnel.     Dimen- 
sions are  such  as  to  make  all  angles  considered  to  be  measured  in 
degrees. 


a     -0° 
o 

a      =   12° 
o 

A 

-    .373 

-      .282 

.364 

B 

+   .0705 

+      .078,7 

+        .0787 

C 

+8.42 

-80.9 

-119.5 

E 

-60.6 

-69.0 

-   69.0 

F 

-1. 237 

-   1.237 

-      1.237 

G 

+1480 

+1480 

+1590 

H 

0 

-133.  5 

+1269 

L 

0 

-133. 5 

-612 

M 

0 

-.0149 

-.0526 

N 

0 

+13.3 
TABLE  A.  1 

+46.8 

AERODYNAMIC  COEFFICIENTS 
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APPENDIX  B 
SYSTEM  CHARACTERISTIC  EQUATION 


B.l    Ideal  System  Characteristic  Equation 

Expansion  of  the  determinant  (3.7)  gives  the  ideal  system 
characteristic  equation  in  symbolic  form 

A    =  1  -  F      FF     -  F  „  F   F     -FF       -FF     -F      F 
+  F2  [  -F8  -  F13  F5  F9  -  Fg  F12  +  Fg  F^  Fg] 

+  F11F13[F5F10-F5F4F9-F7F8F3-F9F7+F8F4F6+F10F5F3] 

+  F11["F10"F12F4F9"F14F8F3"F9F14"F8F4  +  F10F12F3] 

+  FUF1[-F5P10-F12F7F9-F14F8F6+F9F5F14-F8F7+F10F12F6] 
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Upon  substituting  the  aerodynamic  coefficients  evaluated  at 

or    =  24  .  the  C.E.   becomes: 

o  * 

A  =  -  s2[516  s2  +  (50  KL  +  65,  400  K3)s  +  65,  400  K2] 

-  s2(s+l.  237)[-s3+(-69K3  +  .0787  ^  -  .364)  s2  -  (103.5  K2) 

+  (Kg  s2  +  K5  s  +  K4)[-120,  000  s  +  37,  200  KJ 

+  s2  (Kg  s2  +  K5  s  +  K4)[1590  s  +  (-92.  8  K,  +  511.  3)] 


+  (Kg  s2  +  K5  s  +  K4)(K3  s  +  K2)(81,  200  s  +  40,  650) 


B.  2    Non-Ideal  system  characteristic  equation 


(B.3) 


To  obtain  the  characteristic  equation  of  the  coupled  system 
when  non-ideal  sensing  devices  and  non-ideal  control  surface  servos 
are  employed  it  is  useful  to  note  that  these  affect  only  the  "roll  control" 
and  "yaw  control"  equations  of  (3.1). 

The  "roll  control"  equation  is  modified  to 

6(s)  =-G5(KgS2+K5s+K4)G2       (— 2)J(s)  (B.  4) 

s 

The  "yaw  control"  equation  is  modified  to 

S7(s)=  G4[K1(A/3+B7+M6)G1  +  G3(K3s  +  K^]  (B.  5) 

When  these  are  case  in  the  form  employed  in  (3.1)  the 
system  determinant  becomes 
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i  I 


-  sin  a 

-B 

-M 

0 

-E 

-N 

s-F 

-L 

-G 

G2G5X 

2 
s 

0 

1 

s-A  cos  a 

-C  s 

-H  0 

0  0 


(-I^G^A)      (-G3G4Y)  0  (s-I^G^B)      (-KG^M) 

(B.6) 

where  G,  through  Gj.  are  as  defined  in  Chapter  1  Section  1.  2.     This 
constitutes  a  modification  of  the  transfer  functions  defined  in  Chapter 
3  based  on  these  two  equations.     The  following  modified  transfer 
functions  are  therefore  defined  to  replace  the  corresponding  transfer 
functions  of  Chapter  3. 

j.  G9G-X 

j  (.).--%*- .    FUa  (B.7) 

<f>  S 


K,G  G  A 
f (S)  =  jti^G^)    ~=    F12a  (B-  8) 


a  c  y 
f(s)=    (b*K.O,G,B)    S    F13a  <B- 9> 

qj  14     1 


ICGX^M 
P)=     (s-KffiB)    =    F14a  (B'10) 


-K,G  ,G, 
(A/3+B7)^S)s(s-K1G4G1B)  5    F.a.a  (B'  U) 
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With  these  modifications  the  same  block  diagrams,   Figures  3.6  and  3.7, 
can  be  used  and  the  expanded  form  of  (3.  7) ,  see  (B.  l),be  readily  modified. 
For  this  purpose  it  is  convenient  to  regroup  the  terms  of  (B.  1)  as  follows: 

A   "  1  -  F13a(  "P5F3  "  F6  "  F2F5F9  +  F2W 

+  F12a(-F1F6  "  F3  "  F2F9] 

+  FllaF13alF6F10-F5F4F9-F7F8F3-F9F7  +  F8F4F6  +  F10F5F31 

+  Flla!  -F10-F8F4-F1F5F101 

+  FllaF12al  "F4F9  +  F10  F3-FlF7F9+F10FlF6l 

+  FllaF14a[  "F8F3  "  F9  "  F1F8F6  +  F1F9F5' 

"  F1F5  "  F2F8  "  F7F8  (B-12) 


Upon  substitution  for  all  defined  transfer  functions  except  F,, 


through  F.4  there  results 
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Since  the  control  surface  servos  and  the  sensors  are  presumed* 

to  be  second  order  devices,   it  is  possible  to  assess  the  effedts  of  the 

modified  transfer  functions  Fin     through  F,  .     as  follows: 

11a  &        14a 


Define 


G, 


s     +  as  +J? 


(B.14) 


s     +  cs  +  d 


(B.15) 


G, 


2 
s 

+  es  +  f 

h 

2 

s 

+  gs  +  h 

n 

s     +  ms  +  n 


(B.16) 


(B.17) 


(B.18) 


Then  the  modified  transfer  functions  become 


n  X 


(s     +  cs  +  d)      (s     +  ms  +  n)  s 


(B.19) 


KiArz 


(s     +  gs  +  h)      (s     +  as  +  b) 
12a  T^  p  h  b 

s  -  KiB  rs •  -* 


(s     A  gs  +  h)      (s     +  as  +  b) 


(B.  20) 


KxAhb 


2  2 

[s  (s     +  gs  +  h)(s     +  as  +  b)  -  K,Bhb] 


(B.  21) 


♦Telephone  conversation  with  Mr.   R.  L.   Geisberg,CVAC,  Pomona, 
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F 


13a        (s2  +  es  +f)  (s2  +  gs  +  h) 
~~2  ,  v  \    2 


s(s     +  gs  +  h)  (s    +  as  +  b)  -  fr^Bhb 

2  2 

(s     +  gs  +  h)(s     +  as  +  b) 


2 
f  h  (s     +  as  +  b) 


[  s  (s2  +  gs  +  h)(s2  +  as  +  b)  -  K-BhJ?]  (s2  +  es  +  f) 


(B.  22) 


K,M 


h b 


kiiV1  — 2 •     ""5 

„  (s     +  gs  +  h)      (s     +  as  +  b) 

r14a  " 


2  2 

[  s  (s     +  gs  +  h)(s     +  as  +  b)  -  K,Bhb] 

2  2 

(s     +  gs  +  h)(s     +  as  +  b) 


ICMhb; 

— 2 — 2 ; (B«  24) 

[s  (s     +  gs  +  h)(s     +  as  +  ty  -  K,Bhb] 


From  these  it  is  clear  that  the  common  denominator  of  equation 
(B.13)  will  contain  the  product. 

3     2  2  2  2  2 

s   (s   +cs+d)(s   +ms,+n)[s(s   +gs+h)(s   +as+b)-K]Blib]  (s   +es+f)(s-A)(s-F) 

(B.  24) 

3 
Of  these  poles  pf  A  ,   s   (s-A)(s-F)  are  as  in  the  ideal  system.     The 

other  poles  are  due  to  the  factors  of 

2  2  2  2  2 

(s   +cs+d)(s   +ms+n)[s(s   +gs+h)(s   +as+b)  -  K,Bhb]  (s   +es+f)r  (B.  25) 
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APPENDIX  C 


CHU'  S  METHOD 


It  is  the  purpose  of  this  appendix  to  briefly  describe  the  pro- 
cedures by  which  the  "system  determinant"  employed  in  Chapters  2 
and  3  is  derived  from  the  "standard  block  diagram".     A  complete 
development  of  these  standardized  manipulations  as  well  as  a  con- 
siderable body  of  analysis  and  design  theory  for  multiple  loop 

(C  1)  (C   2) 

systems  can  be  found  in  the  works  of  Chu      *     and  Thaler      '      . 

It  is  clear  from  the  manipulation  performed  in  the  early  part 
of  Chapter  3  that  the  principle  of  block  diagram  representation  of  a 
set  of  linear  integro -differential  equations  is  simply  one  of  expressing 
in  LaPlace  notation  the  relations  in  each  equation  as  a  summation 
process  wherein  several  quantities  are  employed  to  form  a  single 
output  quantity  from  a  node.     This  process  is  repeated  for  each 
equation,   and  finally,   interconnections  between  identical  variables  are 
made  to  complete  the  block  diagram.     This  procedure  is  in  common 
use. 

From  the  intermediate  block  diagram  thus  derived,   the  system 

characteristic  equation  can  be  gotten  by  any  one  of  a  number  of 

(C    3) 
methods  including  that  due  to  Mason      *      as  noted  in  Chapters  2  and  3, 

(C   1) 

however,   the  method  developed  by  Chu  provides  additional  facility 

in  performing  analysis  for  each  of  the  several  system  variables  and 
in  design  of  compensation. 
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Figure  C.  1  is  the  block  diagram  of  an  "n  node  system"  (n  equations  in  n 
variables).     The  diagram  permits  of  consideration  of  the  introduction  of 
any  number  of  inputs  (forcing  functions)  to  any  combination  of  nodes.     Any 
block  diagram  derived  by  the  common  process  stated  above  can  be  recast 
into  the  "standardized"  form  by  simply  listing  the  nodes  from  left  to 
right  in  any  sequence  and  performing  the  simple  diagram  manipulations 
necessary  to  accommodate  the  following  conventions: 

a.  Each  nodal  point  may  receive  any  number  of  signals,  but 
delivers  only  one  output. 

b.  Signals  arriving  at  a  node  from  a  node  to  the  left  (feed 
forward)  are  arbitrarily  considered  to  enter  a  node 
additively.     (Sign  adjustment  in  the  feed  forward  transfer 
function  may  be  necessary  to  accommodate  this  convention.  ) 

c.  Signals  arriving  at  a  node  from  a  point  to  the  right  (feedback) 
are  arbitrarily  considered  to  enter  a  node  subtractively. 
(Sign  adjustment  in  the  feedback  transfer  function  may  be 
necessary  to  accommodate  this  convention.  ) 

d.  Signals  to  be  fed  forward  are  derived  at  a  connection 
immediately  to  the  right  of  a  node. 

e.  Signals  to  be  fed  back  are  derived  at  a  connection  immedi- 
ately to  the  left  of  a  node.     (This  and  (d)  frequently  require 
moving  a  "pickoff  point"  where  a  signal  actually  exists  in 
the  basic  block  diagram  to  the  right  or  left  past  a  transfer 
function  block.  ) 

Referring  to  Figure  C.l  and  considering  the  equilibrium  equality 
at  node  a    we  have,  for  the  sign  conventions  noted  above: 

a  =  -  (G  ,  G.    )a  -  (G,G„   )b  - -(G     G,T  )n  (C.l) 

ab    Aa  be    Ba  '  ■  no    Na 

+  G,    I  +  GQ     II  + G      \A 

la  2a  na^' 
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m 
\ 


-v\ 


+  (G     G.T  )n  =  )    (G       )m 

no    Na         i_;        ma 


m=I 


This  equation  can  be  rewritten  as: 


(1  +  GabGAa>a  +  (GbcGBa>  + 


At  node  b  we  have  the  equality: 


b  =  (G  ,  )a  -  (G,     G_,,  )b  -  (G    ,G_  )c  - -  (G     GAT,  )n 

ab  be    Bb  cd    Cb  no    Nb 


(C.la) 


+  G1U I  +  G0,  II  + 
lb  An 


G 


nb 


(C.2) 


which  can  be  rearranged  to  give: 


"  <Gab)a  +  (1  +  Gbc^b  )b   +  (GcdGCb^  +mL^j  '  +  (GnoGNb)n 


■  1    <Gmb)m 


(C.2a) 


m=I 


Extension  of  this  process  to  several  nodes  using  care  to  maintain  the 
sequence  established  in  (C.  la)  and  (C.  2a)  it  can  be  seen  that  the 
sequence  of  nodal  equations  of  equilibrium  can  be  expressed  in  a 
determinant: 


A  = 


^+GaAGAa 


-G 


ab 


-G 


ac 


-G 


an 


G,     G^ 
be    Ba 


-G 


be 


bn 


G    ,G_      ....     G   „G,T 
cd    Ca  no    Na 


1+G,     Gt-,,  G^  jGp,     ....     G^^G-k-t, 

be    Bb  cd    Cb  no    Nb 


1+G    ,G„    ...     G     G-.-T 

cd    Cc  no    Nc 


-G 


en 


1+G     G.T 

no    Nn 


(C.3) 


Inspection  of  (C.  3)  leads  to  the  conclusion  that,  given  the  block  diagram 
in  standard  form  it  is  possible  to  write  down  the  elements  of  the  system 
determinant  by  inspection  of  the  diagram,  using  only  so  much  of  the  n 
node  determinant  pattern  as  applies  to  the  system  at  hand.     The  main 
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diagonal  divides  feedback  and  feed  forward  regions  and  terms  of  the 
main  diagonal  are  the  characteristic  equations  of  the  interior  loops  at 
the  node  indicated  by  the  row  and  column  of  the  main  diagonal  element. 
Expansion  of  the  determinant  in  any  given  case  produces  the  system 

characteristic  equation  without  the  often  confusing  mental  bookkeeping 

(C    3) 

necessary  with  some  other  methods,   including  that  of  Mason      "      . 

By  suitable  choice  of  the  method  of  selecting  the  sequence  of 
nodes  for  drawing  the  "standard  block  diagram"  it  is  possible  (as  in 
Chapter  3)  to  graphically  illustrate  in  the  determinant  the  couplings 
that  exist  within  the  original  equations.     Solution  of  the  equations  for 
desired  input  -  output  relations  follows  Cramer's  Rule,   the  column 
into  which  the  forcing  functions  are  substituted,   representing  the  node 
in  the  standard  block  diagram  at  which  the  desired  variable  exists. 
The  system  being  linear,   superposition  is  applicable. 
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APPENDIX  D 


SCALED  COMPUTER  SETUP  AND  CORRELATION  RESULTS 


D.  1.     Scaled  Layout 

The  REAC  layout  of  the  simulated  system  is  shown  in  Figure 
D.  1.     Note  that  a  time  scale  factor  of  10:1  is  used  (i.  e.  ,   response  is 
slowed  so  that  the  computer  takes  10  seconds  to  respond  as  the  real 
system  does  in  1  second). 

D.  2.     Scaled  Compensator 

The  compensator  V(s)  of  Chapter  3  is  approximated  as  shown 
in  Figure  D.  2.  It  is  limited  in  both  the  inability  to  produce  a  perfect 
derivatiye  and  by  some  approximation  of  time  constants  to  conform  to 
the  RC  values  available. 


Fig.   D.  2    Scaled  Compensator  Layout 
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D.  3    Response  Correlation 

Figures  D.  3  and  D.  4  are  correlation  check  runs  between 
computer  and  analytic  responses  of  the  uncoupled  roll  and  yaw  sub- 
systems. 

Gains  used  are  not  the  final  selected  values  of  equations  (5.1) 
through  (5.6)  but  are  listed  below: 

(D.l) 

(D.  2) 

(D.  3) 

(D.4) 

(D.  5) 

(D.6) 

When  these  values,   analytic  responses  below  are  obtained: 

<f>(t)  =  R[.  464e"i2*76t+  .  465e  ~8'  3Usin(56.  28t-85.  5°)] 

(D.7) 


and 


Kl 

- 

43.8 

K2 

= 

6.18 

K3 

= 

.  357 

K4 

= 

27.9 

K5 

= 

2.  33 

K6 

= 

.0191 

*(t)  =  -P[.  840  +  .  793e"3'  25t  +.  907e"7-  85tsin(16.  63t-106.  4°)] 


(D.  8) 

where 

Pe 


(Aj3+  By  )   I    =  1  ^  (D.  9) 

c  sec 


and 


R=       6T       =ldeg  (D.10) 


In  Figures  D.  3  and  D.  4,   the  continuous  lines  are  computer 
response,   circled  points  are  check  points  from  equations  (D.  7  and  (D.  8). 
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